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The space-time development of a quark-gluon plasma is 
calculated from a Vlasov equation for the distribution func- 
tion of quasiparticles with medium dependent masses. At 
each space-time point the masses are calculated selfconsis- 
tently from a gap equation, whose form is determined by the 
requirement that in thermal equilibrium and for a range of 
temperatures the energy density of the quasi-particle system 
is identical to the one from lattice calculations . The nu- 
merical solutions of the Vlasov equation display confinement. 
Relations to effective theories like that by Friedberg Lee and 
Nambu Jona-Lasinio are established. 

PACS numbers: 25.75.-q, 12.38.Mh, 12.39.Ba 

I. INTRODUCTION 

In view of the possible creation of a quark-gluon- 
plasma (QGP) in high-energy heavy ion reactions at 
the SPS, RHIC and LHC accelerators, descriptions of 
the space-time evolution of the deconfined phase and 
its transition into the hadronic phase are very much 
called for. Ideally, the transport equations should be di- 
rectly derived from the QCD-Lagrangian, a goal which 
is not yet achieved. Instead a number of approximate 
formulations have been proposed. They are either de- 
rived from effective Lagrangians (like those by Friedberg 
Lee |Q] or Nambu Jona-Lasinio (NJL)) or they are 
of Monte-Carlo cascade type |3| with experimental input 
like string-fragmentation functions and/or cross sections. 
Both approaches have their merits. While the ones based 
on effective Lagrangians permit to study certain aspects 
of QCD (e.g. chiral symmetry) and their manifestation 
in the space-time evolution of a highly excited strongly 
interacting system, the phenomenological approaches are 
very helpful in the interpretation of the data. 

The approach presented in this paper is closer in spirit 
to the transport theories based on effective Lagrangians. 
We propose a Vlasov equation for the evolution of par- 
tons, where the medium dependent mass is directly ob- 
tained from QCD, more precisely from the results of lat- 
tice calculations. We concentrate on the transport prop- 
erties of the QGP close to the phase transition, but do 
not describe hadron dynamics. The lattice calculations 
for QCD at finite temperature show a phase transition 
at a critical temperature T c which separates the regime 
of the QGP from the one of the hadrons. Two aspects 
are related to the phase transition: restoration of chiral 
symmetry and confinement. They manifest themselves 
in rapid temperature variations around T c of the chiral 



condensate and the Poliakov loop, respectively. Both as- 
pect are contained in the results (often called " data" ) of 
lattice calculations. 

The basic assumption of our approach is the quasi- 
particle picture, with an effective mass m, which in ther- 
mal equilibrium depends on the temperature T, and on 
space time (a;, t), in the non-equilibrium situation. We 
treat the equilibrium case in Sec. II of our paper and 
show, how m(T) can be obtained from results of the lat- 
tice calculations. Sec. Ill deals with the non-equilibrium 
case, where m{x, t) is calculated from the same gap- 
equation as m(T). We report results of a numerical calcu- 
lation for the expansion of a parton plasma, which shows 
confinement. The relation to other models is given in 
Sec. IV. 



II. THERMODYNAMICS OF THE 
DECONFINEMENT PHASE TRANSITION IN A 
QUASI-PARTICLE APPROACH 

The quasi-particle model is one of the most simple ap- 
proximations to an interacting many body system. A 
system which consists of particles with definite mass and 
interactions among them is replaced by a system of non- 
interacting quasi-particles whose masses m(T) are chosen 
so that the thermodynamics of the original system is best 
approximated. For the case of QCD we make an ansatz 
for the pressure density (thermodynamic potential) for 
the partons Q: 

JVm 1>ma , 1>/ (27r)3 sEftj fi&W 

-V{mi,m 2 , . . .) , (1) 

where the sum runs over the parton species i with degen- 
eracy factor gi and where fi(Ei(p)) are Bose or Fermi dis- 
tribution functions, which de pend on the quasi-particle 
energies Ei(p) — \Jp 2 + mf. V(mx,m2, . ■ ■) describes 
the mean-field contribution to the dispersion relation [Q . 
The potential density V contributes nontrivially to the 
thermodynamic relations whenever the masses are tem- 
perature dependent. 

The masses rrii appearing in the pressure are phe- 
nomenological parameters which have to be chosen so 
that the thermodynamic potential (|l|) has a minimum, 

(|^) =0, i = 1,2,... , (2) 
\drrii J T 

leading to : 
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dV f crp to, ■. , . 

These equations have the form of a gap equation 
for the masses rrii(T) provided the potential density 
V(mi,ni2, ■ ■ ■) is given. Eq. ([|) is equivalent to the 
consistency relations (4|J|] allowing to obtain the energy 
density from the pressure ([!]) , in the form as for the ideal 
gas : 

^p{T) = / T^pEiMEi) + V(T) 

= e kin (m,T)+V(T) , (4) 
where V(T) = V(mi(T), m 2 (T), . . .)• 
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FIG. 1. The energy density of the four-flavor QCD divided 
by T 4 as a function of the scaled temperature T/T c . The 
points are from the lattice data ||. The solid line is the re- 
sult of the quasiparticle model described in the text with the 
chiral symmetry restoration (case (b) and solid line in Fig. 2). 
The short-dashed and the dotted lines are the corresponding 
kinetic energy density and potential energy density respec- 
tively. The dashed line represents the energy density obtained 
as a parameterization of the data points and leading to the 
parton mass increasing at high temperature, (case (a) and 
dashed line in Fig. 2). The dashed-dotted line indicates the 
Stefan-Boltzmann limit. 

We want to model the deconfincment phase transition 
of QCD in a quasi-particle approach and therefore have 
to set a criterion how to determine the quasi-particle 



masses rrii (T) or equivalently the mean-field energy den- 
sity V(T). To date, the most detailed information about 
the QCD phase transition comes from the lattice cal- 
culations which give various thermodynamic functions, 
among them the energy density £i a t{T). We define our 
quasi-particle approach by the requirement 

e QP (T) = e lat (T) . (5) 

Since Eq. ([5]) is only one relation, only one function m(T) 
can be determined. We therefore have to assume that our 
quasi-particle system consists of just one kind of massive 
partons (i = 1 in Eqs. (|l]) to ([!])). We will assume a 
Boltzmann type distribution function f(E, T) — exp _ ^ B . 
The right hand side in Eq. (^) being given, this equation 
is a constraint on the unknown function m(T) in e qp (T). 
We note that the quasi-particle model constrained by Eq. 
(0) describes exactly the energy density of lattice QCD, 
but may reproduce other thermodynamic function only 
approximately |^,|| which is not a surprise since lattice 
QCD is not a quasi-particle gas. The functional form for 
m(T) can be obtained by solving the differential equation 

dm(T) _ ( deiatjT) _ de kin (m,T) \ 
dT ~ V dT dT ) 

I ( de kin (m,T) + dV\ 
\ dm dm J ' 

which is obtained from Eq. (^) by differentiating both 
sides with respect to T and where dV /dm is given by 
the gap equation (||). Eq. (^|) is a first order differential 
equation which determines m(T) for a given energy den- 
sity eiat(T) and for an initial value to (To)- Then V(T) 
can be obtained from V(T) = e lat (T) - e km {m(T),T). 
The lattice data are normalized so that ei a t(T) — > for 
T —> 0. It requires that V(T) — > at low temperatures. 
If to(T) is not obtained from Eq. ^ but is given from some 
other considerations, the gap equation can be integrated 
to obtain the potential V(T) : 

rm(T) f d 3 

V(T) = -g / dm / P f(E(p)) + V(T ) , 

Jm(T ) J (2tiT E[p) 

(7) 

where g = Y.i9i- 

We apply the above methods to the lattice data from 
H shown in Fig. [l|. These data are calculated for 4 
flavors, corrected for finite lattice size and extrapolated 
to massless fermions. We use g = 62.8 corresponding to 
the noninteracting limit of QCD with 4 massless flavors. 
We draw the attention to the fact that the lattice data 
in Fig. |f| for e; a t(T)/T 4 do not approach the Stefan- 
Boltzmann limit for T>T C . This may have two reasons. 

(a) The parton mass m(T) never approaches the chiral 
limit to = 0. Perturbative arguments (whose valid- 
ity is questionable around T c ) suggest that m(T) ~ 
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T for large T. If we attribute the discrepancy be- 
tween Clat (T)/T 4 and the Stefan-Boltzmann limit 
to this reason, one finds the minimal mass m(T) ~ 
2.1T C right above T c and m(T) ~ LIT for large T. 

(b) Chiral restoration requires that m(T) — > above 
the phase transition at least for the fermions. Then 
the deviation from the Stefan-Boltzmann limit has 
to be attributed to a mechanism which is outside 
the scope of the quasi-particle approach. 



T > T c . The mass m a (T) increases linearly with T for 
large temperatures while m&(T) is set to reach the chiral 
limit at high temperatures. Fig. [jjalso shows also the en- 
ergy density corresponding to the two cases. We observe 
that in the case (a) e qp (T) follows the lattice points as it 
should. In particular it has the same high temperature 
limit which is different from the Stefan-Boltzmann limit. 
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FIG. 2. The dashed line represents the temperature depen- 
dence of the quasiparticle mass as obtained from the lattice 
data (case (a)). The short dashed line represents the extrap- 
olation to the low temperature region. It is taken to be con- 
stant, and suffices to effectively confine the partons (an abso- 
lute confinement would require that the mass goes to infinity 
at low temperatures). The solid line represents the assumed 
temperature dependence of the parton mass with zero limit 
at high temperature (case (b)). 

We will explore both possibilities and call the correspond- 
ing masses m a (T) and m&(T) respectively. First we solve 
theEq. (§) to obtain m a (T) for the case (a) , where we set 
m a (T ) = 9.5T C for T = 0.8T C . For the case (b) we take 
m a (T) obtained in (a) for T < T c and extrapolate it at 
T > T c so that m(T) -> at large T (Fig. §). From this 
new functional form mt,(T) we obtain V&(T) and £b(T) 
via Eq. (Q) . Fig. ^ shows the obtained temperature de- 
pendence of the parton mass for the two cases. While the 
two functions m(T) coincide in the region of the confine- 
ment transition (below T c ) they differ dramatically for 
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FIG. 3. The potential energy density V(T) as a function of 
the temperature for the case (a) and (b) of the temperature 
dependence of the mass, dashed and solid line respectively. 



The energy density for the case (b) is different for T > T c 
and approaches the Stefan-Boltzmann limit. Fig. |l| also 
shows the densities for the kinetic energy e/-i n (m(T), T) 
and for the potential energy V(T) for the case (b). Fig. 

shows the potential density functions V(T) for both 
cases. In the case (b), one can identify the large tem- 
perature limit of Vb(T) with the bag constant B, since 
Vb(T = 0) = 0. Indeed for large T (in the deconfined 
phase) Eqs. (|l|) and (^J) take the form of the free mass- 
less gas with a bag constant B = Vb(T = oo). For the 
case (a) one can define the bag constant as the value of 
the potential density V a (T) for the temperature where 
the mass m a {T) is the smallest (slightly above T c ). We 
find 
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£1/4 



1.35T C = 240 MeV 
1.67T C = 300 MeV 



for case (a) 
for case (b) 



(8) 



for a value of T c — 180 MeV corresponding to four- flavor 
QCD. The order of magnitude of the bag constant is cor- 
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rect, but in order to compare it to the usual bag constant 
(B 1 / 4 = 135 - 200 MeV) one should use the lattice re- 
sults for the energy density and the critical temperature 
of the two-flavor QCD as an input for the quasi-particle 
parton model. 

In order to distinguish between the two solutions for 
m(T) we investigate the temperature dependence of the 
chiral condensate in the quasi-particle picture 

(*#) 9P (T) = (W) V ac(T)+g q ,J (0f~/(£) , (9) 

where g qq - counts fermion and anti-fermion degrees of 
freedom and the vacuum part of the chiral condensate 
is given by an expression including a cutoff in momen- 
tum 
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FIG. 4. The chiral condensate as a function of the tempera- 
ture for the case (a) and (b) of the temperature dependence of 
the mass, dashed and solid line respectively. The data points 
are for the lattice data in the four-flavor QCD Q. 

The value of the cutoff A = 430 MeV is fixed to reproduce 
the zero temperature value of the chiral condensate which 
we take twice the usual value 

(**>*«c(0) = 2(250MeV) 4 , (11) 

because we are modeling the four-flavor lattice QCD. 
Fig. |4| shows a comparison between the condensate func- 
tions for the quasi-particle picture and the lattice data. 



This comparison clearly favors case (b), where m(T) — » 
above T c . In what follows, we will mainly work with this 
solution. 

Before we treat the nonequilibrium case we apply the 
gap equation to the case of finite density and calculate 
the density dependence of the parton mass. We introduce 
variables which will be also useful for the discussion of the 
nonequilibrium case. First let us write the gap equation 
in the form 

%-=V(p)^ = -p, (12) 
am am 

where we define the scalar density 
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FIG. 5. The quasiparticle mass as a function of the scaled 
Fermi momentum pf/T c . The solid and dashed lines corre- 
spond to two different temperature dependences of the parton 
mass in Fig. 2. 

It allows us to use the functions V(p) and m(p) instead of 
V(T) and m(T) for the parameterization of the potential 
energy density and the mass respectively. Assuming that 
the energy density depends on temperature T and chemi- 
cal potential p only through p{T, p) we can generalize the 
finite temperature case also to finite density. In principle 
the potential V(p) could depend on other quantities, e.g. 
the baryon density. Any such more general case cannot 
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be discussed using only the lattice data at finite temper- 
ature. As a support for our assumption we note that for 
the NJL model the potential density V depends only on 
the density p (Sec. IV). 

We write the gap equation (0) at finite density 



dV_ 

dm 



d 3 p 



:®(Pf ~ bl) 



(14) 



(2tt) 3 E(p) 

where gf counts the number of fermion degrees of free- 
dom. Note that here we are using the Fermi distribution 
for the fermions at finite density and zero temperature 
and V(m) from the lattice data. Fig. ^ shows the mass 
of the partons as a function of the Fermi momentum pj. 
The behavior is similar as in the finite temperature case. 
For low density the mass increases leading to the effective 
confinement. At high density the mass is proportional to 
Pf for the case (a) and goes to zero in the case (b). 



III. TRANSPORT THEORY OF THE EFFECTIVE 
CONFINING MODEL 

In this section we will discuss the nonequilibrium evo- 
lution of the parton densities. In the semiclassical limit 
the collisonless plasma of quasi-particle partons can be 
described by the Vlasov equation for the phase-space dis- 
tribution function f(x,t,p): 

P 



d t f(x,t,p) 
m(x, t) 



-y x f(x,t,p) 



E(jp,x,t) 
V x m(x,t)V p f(x,t,p) = 



(15) 



f(x,t,p) = -p(x,t) , (16) 



E(p,x,t) 

Eq. ( |l5| ) has to be supplemented by an equation for the 
space-time dependent mass m(x, f) . A sufficient condi- 
tion for the requirement, that the Vlasov equation de- 
scribes the same physics at thermal equilibrium as pre- 
sented in the previous section is that m[x, t) satisfies the 
same gap equation 

dV f d 3 p m(x,t) 

~d^i ~ ~ 9 J (2tt) 3 E(p,x,t) 

where the thermal distribution function f(E) in Eq. (||) 
has been replaced by the nonequilibrium solution of the 
Vlasov equation f(x, t,p) and where V(m) is the same as 
in equilibrium. The solution of the nonequilibrium gap 
equation is then a function m(x, t) of space and time. 
Eq. ( |l6| ) is however not the most general equation which 
reduces to the equilibrium finite temperature gap Eq. 
(||), but one can add to it terms which depend on the 
spatial and time derivatives of m(x,t). We will come 
back to this question in Sec. |Tvj Here we note that the 
choice in Eq. ( pf ) guarantees that the total energy of the 
system 

E(t)=g I d 3 x [ A^E(p,x,t)f{x,t,p) 



d s x . 

J (2tt) 3 

d 3 x/V{p{x,t)) 



is conserved by the evolution according to the Vlasov Eq. 
(|l5|), i.e. the energy of the system is constant. 

We have numerically solved the Vlasov equation to- 
gether with the gap equation for the two cases of the 
functional dependence of the mass m(T) on the temper- 
ature discussed in Sec. ||, but most of the results shown 
relate to the case where the chiral condensate vanishes 
at high temperatures (case (b) in the previous section). 
We have used the test particle method for the solution of 
the Vlasov equation, i.e. we have made the ansatz 



f(x,t,p) = J2s 3 (x-x j (t))S 3 (p- Pj (t)) , (18) 



where the trajectories Xjit) and Pj(t) of the N test par- 
ticles satisfy Hamilton's equations with 



H(x,p) = y/p 2 + m 2 (x, t) 



(19) 



The initial conditions Xj (0) and pj (0) are chosen so that 
given initial densities for matter and momentum are re- 
produced. The initial density is chosen spherically sym- 
metric. Also in the solution of the gap equation the 
spherical symmetry is imposed by angle averaging. 
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FIG. 6. The parton density distribution at different times, 
as obtained from the nonequilibrium evolution of the initial 
fireball (t = 0). 



At time t = the system is described by the density 
profile shown in Fig. ^| and with a momentum distribu- 
tion corresponding to a temperature T = 1.3T C = 180 
MeV. Fig. H shows the density for different times. The 
system expands for t = 2 and 6 fm/c, but then comes 
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back (t = 10 fm/c). The dependence of m(x,t) on the 
radius at different times is shown in Fig. 0. As expected 
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FIG. 8. The component p x of the momentum of a particu- 
lar test-particle taken from the simulation of the time evolu- 
tion of the region of deconfined plasma. 



FIG. 7. The parton mass distribution at different times, as 
obtained from the nonequilibrium evolution of the fireball. 



the mass is small in the interior and increases towards 
the surface. This increase is responsible for the confine- 
ment. Indeed one can show that the equation of motion 
of one particle in the mean-field of the other particles 
conserves the energy of the particle \/p 2 + to 2 . Thus a 
particle cannot leave the region of deconfined plasma, if 
its initial momentum p satisfies: 



2 „ 2 
P < ^va 



(20) 



where m vac is the parton mass in vacuum and to is the 
initial selfconsistent mass of the parton inside the plasma. 
The vacuum mass should be infinite, if the confinement 
is absolute. In our calculation we take the vacuum mass 
equal to ~ 9.4T C (see Fig. ^|), which effectively confines 
the partons at the temperatures discussed here. Thus 
the partons cannot leave the hot fireball if their initial 
momentum is smaller than ~ 9.4T C , which is the case for 
most of the partons at our initial temperature. 

Figs. H and [)] show the time development of the mo- 
mentum p x and the coordinate x of a particular test par- 
ticle. The particle oscillates between the borders of the 
fireball and at the border its momentum is reduced and 
eventually reversed by the action of the force: 



dp to 

—r = Vo-TO 

dt E x 



(21) 



In contrast to the simple particle motion, partons trav- 
eling in bunches may leave into the vacuum region, since 




FIG. 9. The x coordinate of the same test-particle as in 
Fig. I 



the internally created field leads to small masses inside 
the bunch. This possibility is excluded in our method of 



G 



solution, since we require at each time step that the sys- 
tem stays spherically symmetric. This mechanism then 
generates collective vibrations of the surface of the fire- 
ball, when particles are trying to leave the hot region 



ton density inside and thus a small mass of partons. This 
fragmentation mechanism cannot be studied in our spher- 
ically symmetric mean-field theory. It would require a 
description including the fluctuations of the density. 




FIG. 10. The time evolution of the mean square radius 
of the region of deconfined plasma for the case (a) and (b), 
dashed and solid line respectively 
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FIG. 11. The time dependence of the total and potential 
energy of the parton plasma, solid and dashed lines respec- 
tively 



simultaneously. Fig. [T(j shows the time dependence for 
the mean square radius of the fireball : 



(r\t))=g / d 6 x x 



d 3 p 



f(x,t,p) . 



(22) 



The oscillation which can be seen on the figure reflect 
the collective monopole oscillations of the parton density, 
with partons leaving the fireball and reflected back when 
their mass grows. Analogous collective oscillations have 
been observed in the Vlasov evolution of the nucleon in 
the Friedberg-Lee model || . 

As the volume of the fireball oscillates, its potential 
energy will also oscillate growing for large volumes and 
decreasing when the system is compressed. The period of 
the collective oscillations is basically determined by the 
time which a particle needs to travel from one border of 
the fireball to the other (Fig. ||). For massless deconfined 
partons (case (b) in Sec. O) this time is twice the radius 
of the system divided by the speed of light. The total en- 
ergy is of course conserved (see Fig. [ll]) , to the accuracy 
of our numerical solution, and the kinetic (and then also 
the potential energy) has oscillations of the same period 
as the oscillations of the mean square radius in Fig. [l(| 

We mention the possibility of fragmentation of the fire- 
ball into smaller pieces, each of them having large par- 



We state the main result of this section: the Vlasov 
equation together with the gap equation based on lattice 
data shows confinement. For not too high initial tem- 
peratures the fireball stays compact. Of course the sur- 
face of the fireball can oscillate, but the system remains 
bounded. This behavior is in contrast to the free expan- 
sion of the parton gas, where the parton fireball would 
start to expand and cool down faster. In our picture the 
partons are confined, since no hadronization is included 
(see Fig. 



10) 



IV. 



RELATION TO OTHER EFFECTIVE 
MODELS 



The effective models of the QCD using parton degrees 
of freedom are mostly restricted to the fermion sector. 
Thus NJL type models Q have only fermionic degrees of 
freedom with four-fermion interaction. At high temper- 
ature and/or density the quarks are massless or almost 
and at low densities due to the nonzero quark-condensate 
they acquire a finite mass with value around 350 MeV. 
The quarks are not confined in this theory. The NJL gap 
equation for the quark mass can also be written in the 
form of Eq. (H) where dV/ dm is defined by the expression 
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dV_ 

dm 



(m - m ) 
2G 



6N f 



d 3 p 



\ P \<A 



(2tt) 3 E 



(23) 



where G is the four fermion coupling constant, A is the 
infrared cutoff, m is the current quark mass and Nf 
is the number of flavors. Fig. |l2| shows a compari- 
son between the potential extracted from lattice calcu- 
lations and the one from the NJL model. The difference 
is twofold: (i) Shape: While the potentials for Friedberg- 
Lee and NJL model show a clear minimum at the position 
of the constituent mass in the vacuum, the potential from 
our approach drops monotonically to zero, reflecting the 
confinement (infinite vacuum mass), (ii) Magnitude: At 
to = the potential from our approach differs by about a 
factor 5 from the other approaches. This may be partly 
due to the different numbers of degrees of freedom. While 
Friedberg-Lee and NJL refer to a two-flavor quark the- 
ory (no gluons), the lattice calculation is performed for 
Nf = 4 and gluons. The remaining difference may be due 
to quantitative difference between the four-flavor and the 
two-flavor QCD, which goes beyond a simple rescaling of 
the number of degrees of freedom. 
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FIG. 12. The potential energy density V(m) as a function 
of the mass as extracted from lattice QCD (solid line), for the 
two-flavor NJL model (dashed line) and for the Friedberg-Lee 
model (dotted line). We used T c = 140 MeV in defining V{m). 
The large value of the bag constant B = V(m — 0) for the 
potential density from the lattice data can be due to a larger 
number of degrees of freedom in the four-flavor model. 



The Friedberg-Lee model describes fermions coupled to 
a scalar field a(x,t) which plays the role of an effective 
mass and whose dynamics is driven by a potential U(a). 



Of course as in the NJL model the degrees of freedom 
are restricted to the fermions. The Lagrangian of the 
Friedberg-Lee model can be written as: 



— 1 



U(a) , (24) 



where the fermion field operator also carries the fla- 
vor and color indices. The effective quark mass ga can 
include also a contribution from the current quark mass. 

When the expectation value of the a field is infinity, 
the quarks are confined. The same is also effectively true 
if the vacuum expectation value of the a field is very 
large. The gap equation for the quark mass a is given by 
the classical equation of motion for the a field: 



acr + U (a) = -g < > (x, t) 



(25) 



In the case of homogeneous systems Da = (e.g. in 
the mean-field thermodynamics) the Friedberg-Lee gap 
equation is equivalent to the gap equation used in Sec. 
H if Vim) — U{a). In particular the thermodynamical 
energy density discussed in Sec. |l| can be reproduced in 
the Friedberg-Lee model if not for the neglect of the gluon 
degrees of freedom. The difference between our approach 
and the one by Friedberg and Lee rests in the choice of the 
potential. While they assume certain forms, we let V(m) 
be determined by the lattice data. One should note that 
in the homogeneous systems the confining Friedberg-Lee 
model, with the fermion mass m = k(<t) being a function 
of the field er, is also equivalent to our approach after 
a change of variables a — > m. The infinite value of the 
fermion mass means simply in the language of Eq. ( p4| ) , 
that the vacuum gap equation dV/ dm = has a solution 

at 777 = CO. 

The kinetic term for the a makes a difference for the 
case of nonequilibrium or nonhomogeneous systems. In 
the dynamical evolution of a nonequilibrium system in 
the Friedberg-Lee model, the a field is another dynami- 
cal field not related to the local value of the scalar den- 
sity p (|]. The inclusion of the kinetic term for the a 
field however leads to the problem that the value of the 
a field can go negative. The potential U(a) cannot be ex- 
tracted for the negative values of a from the lattice data. 
Moreover in the cases when the mass of the fermion field 
becomes negative, its evolution cannot be described by a 
semiclassical Vlasov equation. 



V. DISCUSSION 

The description of the deconfincment transition in 
heavy ion collisions is a very important theoretical prob- 
lem. A realistic description could allow to define which 
observables are relevant for the observation of the quark- 
gluon plasma formation. A dynamical simulation is 
wished for in order to extract the properties of the plasma 
from the experimental data. Any such approach meets 
difficulties in the description of the confinement and of 
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the hadronization transition when the temperature of the 
deconfined region drops down. In the present work we 
have addressed only a part of this program, namely the 
influence of the confinement on the dynamics of partons. 

The confinement of partons below T c is described in 
a quasi-particle gas model by the increase of the parton 
mass at low energy densities. The temperature depen- 
dence of the parton mass is extracted from the lattice 
data. In the range of temperatures analyzed (0.8T C < 
T < T c ) the mass of the partons increases to values which 
effectively confine the quarks and gluons. The require- 
ment of the restoration of the chiral symmetry at high 
temperatures fixes the temperature mass dependence for 
T > T c . 

The formalism is generalized to the nonequilibrium 
case. The time development of a deconfined fireball stud- 
ied in a Vlasov equation is different from what is usually 
discussed in the literature, in that it shows confinement. 
The expansion of the system is forbidden, but instead we 
observe collective oscillations of the parton plasma. The 
time scale of these oscillations in the collisonless plasma 
is determined by the size of the deconfined region. This 
time scale should be compared to the hadronization time 
in order to determine if the oscillation could develop. 

The increase of the parton mass has implication also 
for the hydrodynamical model of the plasma expansion. 
The slowing down of the hydrodynamic expansion is ob- 
served in simulations of systems with first order phase 
transition Q or in a hydrodynamical calculation with the 
NJL model [[l0| . The use of a confining mass in the hy- 
drodynamical model would stop the expansion. Further 
expansion of the fireball is possible only after hadroniza- 
tion. 

If the hadronization takes place mainly at the surface 
of the deconfined phase then the dynamics of the sys- 
tem could be described by a hydrodynamical model with 
first order phase transition. However, the description 
using the transport equation allows to discuss alterna- 
tive scenarios of the hadronization, e.g. hadronization 
due to parton collisions in the plasma Jllfl with possible 
softening of the spectra of produced mesons. Another 
hadronization mechanism could be the fragmentation of 
the fireball due to a possible instability of the system 
at finite baryon density [H or due to dynamical insta- 
bilities present for energy densities corresponding to a 
mixed phase in the case of a first order phase transition. 
The discussion of the hadronization mechanism and the 
inclusion of the dynamics of mesons remains to be done. 



" on leave from: Institute of Nuclear Physics, Cracow, 
Poland 
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